X-ray Emission from Young Stellar Objects in the e 
Chamaeleontis Group: the Herbig Ae Star HD 104237 and 

Associated Low-Mass Stars 

Paola Testa^'^, David P. Huenemoerder^, Norbert S. Schulz^, Kazunori Ishibashi^'^ 

bo" 

Q : ABSTRACT 

O 
(N 

We present Chandra-RETGS observations of the Herbig Ae star HD 104237 
and the associated young stars comprising lower mass stars, in the 0.15-1.75Mq 
mass range, in their pre-main sequence phase. The brightest X-ray source in 
fH I the association is the central system harboring the Herbig Ae primary, and a K3 

^1 companion. Its X-ray variability indicates modulation possibly on time scales of 

O I the rotation period of the Herbig Ae star, and this would imply that the primary 

'^ ' significantly contributes to the overall emission. The spectrum of the Herbig 

' Ae-|-K3 system shows a soft component significantly more pronounced than in 

other K-type young stars. This soft emission is reminiscent of the unusually soft 
>• ' spectra observed for the single Herbig Ae stars HD 163296 and AB Aur, and 

jT-. ■ therefore we tentatively attribute it to the Herbig Ae of the binary system. The 

O ■ HETGS spectrum shows strong emission lines corresponding to a wide range of 

plasma temperatures. The He-like triplet of Mg xi and Ne ix suggest the presence 



p 

O ' of plasma at densities of about 10^^ cm~'^, possibly indicating accretion related 

oo ■ 

O 



X-ray production mechanism. 

The analysis of the zero-order spectra of the other sources indicates X-ray 
>^ emission characteristics typical of pre-main sequence stars of similar spectral 

^ ■ type, with the exception of the T Tauri HD 104237-D, whose extremely soft 

emission is very similar to the emission of the classical T Tauri star TW Hya, 
and suggests X-ray production by shocked accreting plasma. 
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Introduction 



Nearby young associations provide ideal targets for the study of star and planet for- 
mation, and have sparked growing interest in the past few years. T he early evolution of 



young stars is strongly dependent on stellar mass (e.g., lLarsonlll972l ). For low- mass stars 
(M^ < 2Mq) the accretion timescales are shorter than the evolution time to the zero-age main 
sequence, and these stars cari be observed in their pre-main sequence (PMS) T Tauri phase 
( JD'Antona fc Mazzitellilll994l ). In massive protosta rs (M^ > lOMfr^), however, core hydrogen 
burn ing already starts during the accretion phase (jAppenzellerl Il994l : iBernasconi fc Maeder 
1996). The behavior at the boundaries between high and low mass ranges is yet unclear 
Palla fc Stahk] Il993h . 



Herbig Ae stars (HAe; lHerbig||l960l ) are PMS stars of intermediate mass (~2-10Mq), 
and share several characteristics with T Tauri s tars (TTS), such as IR excess emission, irregu- 
lar ph otometric variability, and disk properties fjWaters fc Waelkenslll998l : [Mannings fc Sargent 
19971 ). HAe stars show evidence for a disk-like geo metry of the circums t ellar material, and 



the disk properties are similar to those of TTS (JMannings fc SargentI 119971 : iGrady et al. 



19991 ). There is some evidence of magnetic accretion in HAe stars, analogous to the magne- 
tospheric accretion scenario generally accepted for classical TTS (CTTS, i.e. TTS that are 
still actively accreting) and are therefore thought to be more m assive analogs of CTTS (e.g.. 



Muzerolle et al.ll2004l : ICradv et al.ll2004 



Guimaraes et al. 



Donati et al 



1997 



20061). While magne tic fields have 



; iHubrig et al.ll2004l ). their origins 



been detected in a few HAe stars (e.g. 

are still unknown; however, a shear magnet ic model does allow ge neration of magnetic fields 

in the early evolutionary stages of A-stars ( Tout fc Pringlelll995l ). 



Systematic X-ray studies have shown that HAe stars are moderately bright X-ray sources 
toamiani et"aLlll994l : Izinnecker &: Preibischiri994l : iHamaguchi et al.lbood : Istelzer et aDl2006h . 
Since main sequence A-stars are not known to be luminous X-ray sources due to lack of strong 
enough winds or magnetic-dynamo driven coronae, the findings may imply that the physical 
characteristics of young A-stars stars differ from those of main sequence A-stars. 

HAe X-ray emission mechanisms remain unknown. Early studies based on ROSAT and 
Einste in observations proposed the standard wind shoc k model as for early-type massive 
stars fJDamiani et al.l Il994l : IZinnecker fc Preibischl Il994l ). However, typical Lx/I/boi ratios 
found for HAe s tars are significantly higher than the empirical luminosity ratio of 10~^ for 
early type stars (jSkinner et al.ll2004l ). casting some doubt for this paradigm. Also, the high 
X-ray te mperature (> 10^ K) component found for a large fraction o f HAe X-ray sources 



[see e.g., ISkinner et al.ll2004j : IHamaguchi et al.ll2005l : IStelzer et al.ll2006l ) cannot be produced 
through this mechanism. Others works attribute the emission to magnetic activity (either 
coronal activity as in late-type stars, or magnetic activity due to star-disk interaction). 
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or to unresolved late-type companions (JHamaguclii et al.l l2005l : ISkinner fc Yamauchil Il996l : 
Skinner et allbood : Istelzer et aPboOGJ ). 



Another emission mechanism possibly at work in HAe stars is X-ray emission from 
shocks in accreting plasma, proposed as the pro duction mechanism for soft X-rays in several 
CTTS (TW Hya, BP Tau, V4046 Sgr, MP Mus 



Kastner et al. 



Schmitt et all 1200,4 lOiinther et al.ll200fil : lArgiroffi et al.ll2007l ). For these CTTS, Chandra 



2002; Stelzer fc Schmitt 2004 



and XMM-Newton high resolution spectra have revealed peculiar characteristics with respect 
to all other observed stellar spectra, in particular showing (i) soft excess — prominent cool 
component at temperatures of ~ 3 x 10^ K, TW Hya being an extreme case, where this 
soft component largely dominates the X-ray spectrum; (ii) unusually high density 



in. 



> 



10^2 cm-3 



for the cool (T ~ 2 — 4 x 10^ K) plasma, as diagnosed through the analysis 
of the He-lik e triplets line ratios ; (iii) peculiar abundances — extremely high Ne, low metal 
abundances (JKastner et al.l |2002| ; IStelzer &: SchmittI |2004| ; iDrake et al.l l2005l ) . Analogously, 
Chandra imaging spectra of the HAe star HD 163296 revealed a very soft spectrurn (~ 
0.5 keV) possibly sugges ting accretion domi n ated X-ray emission (ISwartz et al.ll2005l ). as 
for the CTTS TW Hya JXastner et al.ll2002l ): lOradv et all (|2007l ) also find extremely soft 
X-ray emission (~ 0.25 keV) from the Herbig Ae star HD 169142. Recently, XMM-TVewton 



obser vations of the HAe star AB Aur have shown a very similar soft spectrum (JTelleschi et al. 



20071 ). In this case, the high resolution spectrum provides additional diagnostics: the high 
forbidden/intercombination ratios {f /i) of the He-like triplets, especially O vil, indicate the 
absence of high density and high UV field, expected if X-rays are produced in the accretion 
shocks, and rule out accretion-related X-ray emission for AB Aur. For AB Aur, which has 
spectral type AO, the photospheric UV emission is strong enough to alter the O He- like triplet 
up to several stellar radii from the surface; the large measured O Vll f /i line ratio, constrains 
the emission to originate high above the photosphere (> 3R^) and therefore it implies that 
also coronal emission is unlikely for AB Aur. A plausible alternati ve for this source is 
X-ray emission from magnetically confined winds, as in the scenario of iBabel fc Montmerle 
( 119971 ). It is important to keep in mind that some of the different X-ray emission mechanisms 
might not be m utually exclusive and instead they might coexist, as seems to be the case 
for some CTTS dSchmitt et aPbood ; Iciinther et aPbood ; lArgiroffi et aDl2007h where both 
magnetically confined coronal plasma and accretion streams likely contributes to the overall 
X-ray emission. 



HD 104237 is one of the nearest known {d = 116 pc; iPerryman et al.lll997l ) and well stud- 



ied HAe stars, and it is also a strong X-ray source (Lx ~ 2 x 10^*^ erg s~^). Chandra observa- 
tions have spatially resolved the X-ray emission from HD 104237, showing that the brightest 
X-ray source is coincident with the HAe source, and identifying 4 previously unknown nearby 



sources, likely TTS forming a young stellar group associated with HD 104237 fJFeigelson et al. 
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20031 . hereafter FLG03; see §2]). Chandra (FLG03) and XMM- Newton flSkinner et al.ll2004 ) 
imaging spectra of this HAe star have shown the presence of hot plasma (> 2 keV), taken 
as evidence of magnetic activity. HD 104237 with its proximity, high Lx, and low line-of- 
sight extinction is a unique target to study the X-ray emission mechanism in HAe stars, 
through high resolution spectroscopy. High resolution spectra provide plasma diagnostics 
which are the only means for detailed study of the physical conditions of the emitting mate- 
rial. Cool, dense plasma (logT[K] ~ 6.5, logne[cm~'^] > 12) is now thought to be a signature 
of accretion, while variable hot emission is indicative of magnetic activity. 

In this paper we present Chandra-UETGS observations of HD 104237. The relatively 
long exposure of ~ 145 ks provides detailed information on the X-ray variability of the HAe 
star, as well as of other stars in the field, and it yielded a well exposed hetgs spectrum of 
the HAe star providing detailed diagnostics for the emitting plasma. 

We present the characteristics of the targets in §21 The observations and our techniques 
of line flux measurement and spectral analysis are briefly described in ^ The results are 
presented in §11 We discuss the results and draw our conclusions in §3 and El 



2. HD 104237: the Herbig Ae star and the T Tauri association 

HD 104237 is the optically brightest HAe star, and it is a lso a moderately s trong X- 

ASCA 



ray source with logLv ^ 30.5 erg/s, observed with ROSAT (lAlcala et al. 



(jSkinner fc Yamauchilll996l ). Chandra (FLG03), and XMM- Newton (jSkinner et al. 



11995 



200J). 



Chandra imaging observations, thanks to the unprecedented sub-arcsecond spatial res- 
olution, have revealed the presence of four other X-ray emitting young stellar objects as- 
sociated with HD 104237, labeled B-E, A being the HAe star (Figured see also FLG03). 
Grady et al.l (120041 ) (hereafter G04) have conducted intensive high angular resolution mul- 
tiwavelength imaging (optical, NIR, mid-IR), and spatially resolved spectroscopy (optical, 
UV, FUV) to probe the characteristics of A, and the associated star-forming environment. 
The comprehensive study by G04, has provided detailed information on the HAe star and 
constraints on possible close companions. 

The stellar parameters of the members of this young (~ 2-5 Myr) stellar group are 
listed in Table [H IR excess. Ha emission, and Li absorption indicate a T Tauri nature for 
the components B, D, and E (FLG03; G04). FLG03 hypothesized that C is a late M-type 
TTS or a brown dwarf. 



HD 104237- A: HD 104237-A has been identifled as an Herbig Ae star bv lHu et all (119891 ). 
based on its optical spectrum and IR excess. It was initially assigned a spectral type AOpe by 
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Fig. 1. — Color coded Chandra image of the HD 104237 field from our ~ 145 ks hetgs 
observation. The X-ray sources of the young stellar group are labeled following Feigelson 
et al. (FLG03). The separation between sources is hsted in Tabled} all sources are within 
~ 20". 



Walker fc Wolstencroftl (Il988l ). then A4 ve by lHu et al.l (Il99ll ). and recently G04 have revised 
its classification as A7.5ve-A8ve. The HAe star presents evidence of activity-as indicated 
by line emissi on, UV excess, outflows- and is still accreting at an estimated rate of ~ 10"^ 
Mq/ji ('G04. iGarcia Lopez et al.l 120061) . UV high resolution spectra show P Cygni profiles 
for several lines (JBlondel fc Djiell2006l ). G04 also found a bipolar microjet and Herbig-Haro 
knots; from the analysis of proper motion of a knot they derived the inclination of the system 
to be 18°!}^. The inclination may be a meaningful parameter for the detection of accretion 
signatures, assuming the geometry of magnetospheric accretion models; for instance, the 
CTTS TW Hydrae and BP Tau, whose soft X-ray emission appears accretion dominated, 
have both nearly pole-on aspects of 10° and 30° respectively. The observations suggest a 
magnetospheric accretion scenario for HD 104237-A si milar to CTTS (G 04). A magnetic 
field of this HAe has been marginally detected (~ 50 G; iDonati et al.lll997l ). If the geometry 
of the accretion is similar in HD 104237-A we will likely observe almost along the accretion 
stream. The HD 104237-A disk extends at most to 0.6" (~ 70 AU). 

ifST observations put tight constraints on the distance of any unknown close companion 
to A (other than the resolved X-ray source B at ~ 1.4" separation; G04), excluding the pres- 
ence of a companion with distance r to the A star 0.05" < r <1", imp l ying that a possible 
companion would have to be closer than ~ 5 AU to A. iBohm et al.l (120041 ) through high 
resolution spectral observations (aimed at the study of 5 Scuti pulsation observed for this 
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HAe star) revealed the presence of a spectroscopic binary companion to the HAe star, and 
determined the orbital parameters. They find a mass ratio Mprimary/Msccondary = 1.29 ±0.02, 
and an orbital period of about 20 days, implying an average separation of ~ 0.15 AU. They 
assign a spectral type K3 to this spectroscopic companion which is responsible for the Li- 
line traces in the spectrum of the HAe star revealed by FLG03 and which they attributed to 
contamination from the adjacent B component. Th e parameters of the HAe star and the K3 



close companion are included in Table [2l Recently, iTatulli et al.l (120071 ) using near-IR inter- 
ferometric observations of this HAe star modeled the spatial distribution of the circumstellar 
material and outfiows on AU scales; their observations are compatible with outflowing wind 
launched in the vicinity of the dusty disk inner rim (~ 0.2 — 0.5 AU). Considering the very 
small distance between the spectroscopic binary components, the circumstellar disk of the 
HAe star is likely a circumbinary disk, and this could in principle affect the magnetospheric 
accretion scenario and other phenomena occurring on small scales close to the stars. 

The presence of this later type companion, unresolved in the X-ray observations, sig- 
nificantly changes the perspective for interpretation of the X- ray emission of HP 104237-A, 



as K-type PMS stars are known strong X-ray emitters (e.g., iGetman et al.ll2005l ). We will 
thoroughly discuss this issue in §3 

HD 104237-B: The proximity of HD 104237-B to the HAe star hampers the precise 
determination of its stellar parameters, as thoroughly discussed by G04. FLG03 tentatively 
assigned a K spectral type to B due to the presence of Li-line features in the optical spectrum 
of HD 104237-A; however, these features are actually produced by the the K-type close 
companion subsequently discovered as discussed above. G04 finds that HD 104237-B is 
likely a M-type TTS on the basis of its infrared excess, and of its PSF in HST/STIS images, 
compatible with M3-M4 spectral type. 

HD 104237-C: FLG03 has tentatively classified HD 104237-C as a brown dwarf candidate. 

HD 104237-D: The Ha equivalent width, close to the classical TTS limit, the Li features, 
and the X-ray luminosity indicate the T Tauri nature of HD 104237-D (G04). 

HD 104237-E: This star has a late spectral type, broad Ha emission, and photometric 
variability, which are typical of accreting CTTS (G04). This source was also variable in 
X-rays and showed a remarkable change in emission level by a factor ~ 8 between the two 
short ACIS-I observations of FLG03, who also found it to be much more absorbed than the 
other stars in the association. This suggests that the obscuration is due to local circumstellar 
material. 
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3. Observations and Analysis 



We observed the Herbig Ae star HD 104237 and the group of young stellar objects 
associated to it with the Chandra High Energy Transmission Grating Spectrometer (see 



Canizares et al.l |2005| for a description of the instrumentation) for a total exposure time of 



145 ks. The details of the observations are presented in Table Ei 

We also re-analyzed the ACIS imaging-mode observations of FLG03 (ObslD 2404 and 
3428). 



The data used here have been reprocessed using standard CIAO (JFruscione et al.ll2006l ) 
v3.4 tools. Effective areas were calculated using standard CIAO procedures, which include 
an appropriate observation-specific correction for the time-dependent ACiS contamination 
layer. 

Spectral analysis of grating and zeroth order s pectra were done with t he Interactive 
Spectral Interpretation System (ISIqj) version 1.4.2 ( Houck fc Denicolall2000l ). 



In Figure [T] we present a color coded image of the Chandra-UETGS observations where 
photons energies in the range 2.7 — 7.0 keV correspond to blue, 1.2 — 2.7 keV to green, and 
0.5 — 1.2 keV to red. The X-ray sources are labeled following the classification of FLG03. 

In the following section we present the analysis of the X-ray emission of all sources of 
this young stellar association. We analyzed the low resolution spectra and variability for all 
group members, and the high resolution spectrum of source A which has enough counts for 
a meaningful analysis of the dispersed photons. 



4. Results 

In the total ~ 145 ks hetg observations we detected all sources studied in FLG03. 
Source A is the strongest X-ray source but even for this source the pileup in the zeroth order 
is essentially negligible (up to a few %). Source E and D have comparable count rates about 
an order of magnitude smaller than for source A. Source B, which is at an angular distance 
of > 1" (FLG03, G04), is apparent as an elongation of source A in the SW direction; even 
though there is a small contamination of source B in the spectrum of source A its contribution 
to the A-l-B spectrum is lower than 10%. Finally, source C is detected at a 3a level only for 
a portion of the observing time. 



ISIS is available at http://space.init.edu/cxc/isis/ 



The HETG observations presented here provide significantly better statistics than the 
previous ACIS-I observations analyzed by FLG03. For source A, we obtained more than 3700 
integrated zero-order counts vs. 640 of the short ACIS observations. For the sources B, C, D 
and E we obtain 152, 54, 251, and 543 zero-order counts which are a factor 2.5 to 9 times 
higher than the counts of the ACIS-I observations which were 44, 6, 70, 225, respectively. 

In order to extract the zero-order spectrum of source B we needed to take into account 
the contamination of source A. This contamination was estimated by extracting the counts 
of source A in an annular region 2-4 pixels from the center, and excluding a sector containing 
source B (between -0.64 and 0.24 rad, where 0.0 is west and positive angles are in the north 
direction). These counts were then scaled to the region of extraction of source B and used 
as a background for this source. 



4.1. Variability 

All sources in the HD 104237 field appear to be variable. In Figure [2] we show the 
lightcurves obtained from zero-order photons for sources B, C, D, E. In order to show how 
significant the detection of source C is in the different segments of observation we also derived 
a lightcurve for the background, using an extraction region, close to source C, and with radius 
4 times larger than the default value used for the extraction of the spectra of all sources, and 
then scaling the rate to the source region areas. All sources present significant variability 
over time scales of days to months. 

In Figure [3] we plot the lightcurves of source A as obtained from the dispersed photons. 
In the upper panels we plot the lightcurve of all dispersed photons (black curves), and also 
the lightcurves in two spectral bands in order to show possible changes in spectral hardness: 
the hard band corresponding to the 1.5 — 8 keV energy range (blue curves), and the soft 
band corresponding to the 0.5 — 1.5 keV range (red curves). The lower panels show the 
corresponding hardness ratio (defined as [hard-soft]/ [hard -fsoft]). These plots show that the 
source presents continuous variability (within a factor ~ 2 from an average value) on time 
scales from hours to days. However the observed variability does not show obvious similarities 
with very dynamic events such as flares, which are typical of active stellar coronae and are 
generally characterized by significant increases in hardness ratio and time scales of hours up 
to days for the most extreme events observed in very active young stars. 

We searched for possible periodicity in the lightcurves of source A, in particular exploring 
possible periodic variation on the characte ristic time scales o f the system: the orbital period 



of the binary system (Porb = 19.859 days; iBohm et al.l 120041 ). and the rotation period. The 
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Fig. 2. — Zero order lightcurves of sources B, C, D, E, and for the background: we plot count 
rates with error bars for the four observations as a function of time; the last observation, 
carried out two months apart the first one is shifted by -50 days. The background rate was 
derived from a circular region close to source C, with a radius 4 times larger than the default 
extraction radius used for all the sources. The counts for source B are derived as described 
in the text with a procedure we devised to take into account and subtract the contamination 
of source A. 
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Fig. 3. — Lightcurves of first order dispersed pliotons of source A for each of the four 
observations, using a binsize of 2 ks. The upper panel of each plot shows the hghtcurves 
of all the first order dispersed photons (solid line with error bars), and the lightcurves in 
two energy bands: a hard band (1.5 — 8 keV; dotted line) and a soft band (0.5 — 1.5 keV; 
dashed line). The lower panels show the corresponding hardness ratio defined as [hard- 
soft] / [hard+soft] . 



rotation period of this Herbig Ae star is however not well d etermined, th e refore we used 
different estimates. Using the vsini = 12 km s~^ obtained by iDonati et al.l (119971 ) (derived 
from several photospheric spectral lines), the inclination i = 18° estimated by G04, and a 
stellar radius R^ = 2.5Rq ( iBohm et al.ll2006l ) we obtain a rotation period Prot = 3.26 days; 
however, con sidering the uncerta inties on all these parametercl, the rotation period is poorly 
determined. iBohm et al.l (120061 ) present a possible indication of rotational modulation by 
studying modulation of lines forming close to the stellar photosphere: they find modulation 
in the Ha line that might be due to rotation and it would yield Prot = 100 ± 5 hrs, i.e. 



■^For instance, iBlondel fc Djid ( 20061 ) assume very different values for both f sini (150 km s ^), and i 
(53°), which they infer from optical lines, and the redshifted component of the Lya emission respectively. 
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-Prot = 4.17 ± 0.21 days. Therefore we performed a period analysis in this range of periods 
and we find a best fit period of 95 hrs. 




Fig. 4. — Lightcurves of dispersed photons of source A (binsize of 2 ks) phase-fol ded assuming 
three different periods: top - P = 75.6 h obtained assumin g v sin i = 12 km s~^ fJDonati et al. 
1997h . i = 18° (G04), and R^ = 2.5Rq jB5hm et ahlbooeh: rmddle - P = 95 h (within range 



of possible rotational modulation observed in Ha by iBohm et al.l 120061 : see text); bottom 
P = 19.8 59 days = P^rh o f the binary system of the Herbig Ae star HD 104237 + K3 
companion ( Bohm et al.ll2004l ): phase corresponds to the periastron passage. Filled symbols 
indicate the zero-order count rate predicted from modeling the ACIS-I observations (see text 
and Table Hj). The dotted line represent a fit to the data with a simple sinusoidal curve plus 
a constant; while this function does not represent a specific physical model, the constant 
function can be interpreted as the emission level of the quiescent (i.e. more homogeneously 
distributed) corona, while the rotational modulation would be likely caused by the appearing 
and disappearing from view of compact regions of more intense X-ray emission. 

Figure m shows the phase folded lightcurves for three different periods. All plots indicate 
that our X-ray observations do not provide a large phase coverage even for the shortest period, 
therefore making difficult to identify possible periodic modulation of the X-ray emission. 
Nevertheless, the observed lightcurve seem to be compatible with modulation on time scale 
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comparable with the rotation period as estimated by iBohm et al.l (120061 ): if true this would 
imply that the Herbig Ae star is likely a significant contributor to the overall X-ray emission 
of the binary system. The variability does not show correlation with the orbital period of 
the HD 104237-A system of 19.859 days (bottom panel of Fig. S]). 

In Figure H] we also superimpose the count rates estimated for the two early ACIS-I 
observations studied by FLG03 (ObsID 2404 and 3428). FLG03 in their analysis assumed 
negligible A^h and did not account for pileup. In order to take into account these effects and 
estimate the expect ed count rat e for these two observations, we re-analyzed the data using 



the pileup model of (lDavisll200ll ) as implemented in ISIS. This model accounts for coincident 
detection of multiple photons and has the advantage of using all the counts in the source 
region, not just those in an outer, un-piled annulus. We found a two-temperature model was 
necessary, and that the pileup fraction was about 20%. In fitting the spectra we fixed the 
abundances to the values derived from the hetgs analysis, and the Ah to 2 x 10^^ cm~^ (see 
detailed discussion in §4.2p . We note that our analysis of the ACIS-I spectra of HD 104237-A 
yields for both temperature components lower values with respect to those found by FLG03; 
also, our derived luminosity values are larger than the values found by FLG03 (by about 20% 
and 50% respectively in the two observations). These discrepancies are due to the inclusion 
of interstellar absorption and pileup effects. 



4.2. Low resolution spectra 

The color coded image (Figure [1]) shows that source E is characterized by a hard spec- 
trum, source D is on the contrary extremely soft, and source A is an intermediate case. Here 
we present the low resolution spectra of source A, B, C, D, and E and the results of the fits 
to these spectra with isothermal or two-temperature plasma models. 

For source A, the higher signal to noise allow us to study separately the low resolution 
spectra of each observation to highlight possible changes in the spectral characteristics with 
changes in the X-ray emission level. For all other sources due to the poorer statistics we 
combined the spectra of all four observations in order to improve the signal-to-noise ratio 
and better constrain the average spectral parameters. 

We extracted the low resolution spectra of source A for each observation, and performed 
a fit with a thermal model of plasma in coUisional ionization equilibrium with variable 
absorption. A single temperature model, even with variable abundances, does not provide 
a good fit to the spectrum. Therefore we use a 2-T component model with single Ah- The 
fits to the spectra obtained for the separate observations do not show significant changes in 
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the spectral parameters within the uncertainties, therefore we summed up the spectra for all 
four observations to better constrain the parameters of the model. 

As thoroughly discussed in §4.3.11 the high-resolution spectrum provides additional 
strong constraints to the fit of the X-ray spectrum. Here we use the abundances found from 
the analysis of the dispersed spectra and we fit the low-resolutions spectrum with a 2-T 
model with abundances fixed to those values. This 2-T model, whose parameters are listed 
in Table 0, provides a good fit to the integrated zero-order spectrum of source A, and shown 
in Figure [5l We also fitted the model with a similar 2-T model but with variable abundances 
(model II in Table [5]), however the comparison with model I indicates that this model does 
not provide a significantly better fit: when abundances are let free to vary their best fit 
values are very similar to the values found from the analysis of the dispersed spectra which 
provide more robust results because single lines are resolved and used for the determination 
of the element abundances (see §4.3.11) . 

Furthermore we explored two more different models: a 2-T model with fixed abundances 
as in I but different N^ for each temperature component (model III), and 3-T model with 
fixed abundances as in model I (model IV). The former model explores the scenario of 
two temperature components corresponding to different X-ray emitting plasma regions with 
different N^, as likely for ex ample if the two components are emitted by the two binary stars. 



Skinner Sz Yamauchil (119961 ) already investigated this possibility when analyzing the ASCA 
data and found a hot more absorbed component (T ~ 20 MK, A^h ~ 1-7 x 10^^ cm~'^) and a 
cool slightly less absorbed component (T ~ 2.5 MK, Nn ~ 6 x 10^^ cm"'^). They tentatively 
interpreted the hot absorbed component as produced by an embedded TTS companion, and 
the cool emission as the HAe star emission; the uncertainties however were too large to 
provide any real constraint to this model. From our data we find that the absorption of the 
hot component might be actually lower than the absorption of the cool component, but they 
are not significantly different within the uncertainties. The 3-T temperature best fit model 
has a cool component with parameters essentially identical to the 2-T models, while the hot 
component is redistributed in two hot components at 20 and 80 MK respectively; the only 
parameter slightly different is the lower absorption Ah- Neither of these two model provides 
a statistically better fit to the observed spectrum than the 2-T model with fixed abundances, 
therefore we consider model I (Table E]) the best fit to the low-resolution spectrum. 

Comparing these results with the finding of analysis of previous X-ray observations 
of HD 104237 we note that the overall characteristics (T, Lx) of its X-ray emission seem 
remarkably constant for more than a decade, with Lx ~ 2 x lO^'^ erg s~^ and at least 



two ter nperature components needed to expla in its spectrum (ISkinner Sz Yarnauchi 



FLG03. ISkinner et al.ll2004l : IStelzer et al.ll2006h . Both FLG03 and IStelzer et al.l fl2006f ) find 



1996 
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Fig. 5. — Zero order spectrum of HD 104237-A (data points) with best fit model super- 
imposed (model I of Table [Sj solid line). The best fit model is a two-temperature plasma 
model with Nu= 1.7 x lO^i cm'^, Ti = 4.7 x 10^ K (0.4 keV), T2 = 2.7 x 10^ K (2.3 keV). 
The X-ray luminosity from this model is 2.9 x 10'^'' erg s~^ (corrected for absorption) in the 
0.5-7 keV range. 
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from the analysis of Chandra- AClS data slightly higher temperatures for both T components; 
however this can be explained considering that they assume essent i ally n egligible N^. Our 
model parameters compare well with the findings of ISkinner et al.l (120041 ) based on XMM- 
Newton data, t he on ly significant difference being in the derived element abundances which 
Skinner et al.l (120041 ) find to be significantly above solar for all elements except S, and Fe, 
about solacl. 

The fits to the low resolution spectra of HD 104237-A all yield A^h values compat- 
ible with the visual extinction derived for this star: the T V h val ue spans the range 0.1- 
0.2x10^^ cm~^, which, using the conversion of iGorensteinI (119751 ). corresponds to Ay in 
the range 0.45-0.9, in good agreement with values found fr om optical studies (e.g., 0.3 of 
van den Ancker et ah 19981. and .8 of iMalfait et al.lll998l ). This is also consistent with 
previous findings of ISkinner et al.l (120041 ). 



The zero-order spectra of the other four sources are shown in Figure [6] and [71 and the 
parameters of the fits with corresponding uncertainties are listed in Table [61 

The left panel of Figure [HI shows the zero-order spectrum of source B with the corre- 
sponding best fit model. We find that a model with a single temperature at ~ 1 keV and 
zero absorption reproduces well enough the observed spectrum. The parameters of the X- 
ray spectrum are rather typical of similar M-type TTS (see discussion in ^. However the 
relatively low number of counts (~ 150) obtained for this source does not allow to constrain 
the parameters very well, as demonstrated in the confidence contours in the Ah-T parameter 
space shown in the panel below the spectrum plot in Figure [61 

Source C is a very weak source for which very scarce information is available. On the 
basis of the X-ray observation and the (lack of) optical data, FLG03 proposed a possible 
identification as a brown dwarf; the association with the young stellar group of HD 104237 
is uncertain. In our 145 ks observation we gathered about 54 net counts for source C 
(subtracting a background estimated on a nearby extraction region); the spectrum and the 
best fit model are shown in Figure [61 The poor statistics do not provide stringent constraints 
on the spectral properties of the X-ray emission of this source. Within the uncertainties, we 
find parameters compatible with the emission by a brown dwarf at the distance of HD 104237 
confirming the findings of FLG03. 

For source D, and E, which are identified as T Tauri stars associated to HD 104237, we 



^They find Ne, Mg, Si, Ca, to be 2.3, 2, 2, 9.9 times solar respectively for their 2-T model, and 2.2, 1.7, 
1.6, 6.6 times solar for their 3-T model. As discussed in i i4.3.1l the analysis of the high resolution spectra 
allows to determine the element abundances accurately, ruling out such large values. 
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Fig. 6. — Top: Zero-order spectrum of HD 104237-B {left) and C{right), rebinned to a 
minimum 5 counts per bin, with corresponding best fit model. The parameters of the best 
fit models are listed in Table [61 Bottom: Confidence contours in the N^-T space (68%, 90%, 
98%) for source B {left) and C {right). 
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have obtained about 250 and 540 zero-order counts respectively, over the 145 ks observation. 
As already apparent in the color coded image (Fig. [1]), source E is heavily absorbed and only 
the hard photons are detected. The spectral fit indicates in fact a large N^ of 3.4 x 10^^ cm"^, 
much larger than the interstellar absorption and that can be therefore ascribed to local 
absorption from circumstellar material. This source is identified as a K3 classical T Tauri 
star, and the infrared excess detected by G04 provide further evidence of presence of large 
amounts of circumstellar material. Source D is a lower mass T Tauri star, of spectral-type 
M, border-line classical T Tauri star as indicated by the Ha equivalent width very similar to 
source E (G04). At variance with respect to source E, infrared measurement do not indicate 
large infrared excess (G04; however the measurements are characterized by poorer statistics 
for this source and therefore the limits are not very stringent). Its X-ray spectrum is very 
soft for its emission level (its Lx= 10^^ erg s~^ is rather large for M-type TTS; see §3 for a 
more detailed discussion), and it is fit well with an isothermal model with T = 3 x 10^ K 
(0.24 keV). These spectral parameters for source D, and E, are well determined within the 
statistical uncertainties as clear from the confidence contours in the A^h-T space shown in 
Figure [7| (bottom panels) . 



4.3. HETGS spectrum of HD 104237-A 

A goal of our observing program was to obtain a high resolution X-ray spectrum of the 
Herbig Ae star with high signal-to-noise ratio, in order to investigate the X-ray production 
mechanisms in intermediate mass pre-main sequence stars. The HETGS spectrum of sources 
A and B cannot be resolved spatially in the dispersed spectra. However, our zero-order 
analysis ( §4.21) indicates that component B's contribution is < 10% that of A. Henceforth 
we will refer to the spectrum as if from A alone, but it should be kept in mind that there 
is some level of contamination from B. Other sources in the field are distant enough that 
order-sorting results in unique determination in MEG-HEG crossed-order regions. In the 
145 ks HETGS observation we obtained slightly more than 5000 total counts in the HEG and 
MEG -1,+1 orders. 

In Figure[8]we show the high resolution spectrum of HD 104237-A, obtained by summing 
HEG and MEG spectra (-1 and -|-1 orders). The spectrum is dominated by line emission of 
H-like and He-like ions of Si, Mg, Ne, by Fexvn transitions around 15A and 17A, and by 
transitions of H-like oxygen. Also, some continuum emission is present as well as transitions 
of highly ionized Fe around ~ 11 A both indicating the presence of significant amounts of 
hot (T > 2 X 10^ K) plasma. 

The measured line fluxes together with the statistical errors, are listed in Table [71 
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Fig. 7. — Top: Zero order spectrum of HD 104237-D (left) and E (right) with best fit model. 
Tlie parameters of the best fit models are listed in Table [6l Bottom: Confidence contours in 
the A^H-T space (68%, 90%, 98%) for D (left) and E {right). 
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Table 1. Characteristics of HD 104237 members. 



Object Spec. Type ^ EW(Ha) 



offset 



comments 



Ad 


A7.5Ve-A8Ve 


-24. 


B 


M3-4 




C 






D 


M2-3 


-9 


E 


K3IVe 


-9 



HAe, IR and UV excess, jets 
f .365"±0.019 IR excess 

5.275"±0.033 BD (FLG03) or M6V (G04) 
10". 72 ± 0.05 Ha near CTTS limit 

14". 88 ± 0.05 Ha near CTTS limit, IR excess 



Ranges of spectral types from iFeigelson et al.l (|2003r ): iGradv et al.l (|2004| ): iLuhman 



daooj). 

^Equivalent width, in A, of Ha emission from lAcke et al.l ( 2005 ): lGradv et al.l (2004) 



'^A has a K3 spectroscopic companion within ^ 0.15 AU (JBohm et al.l 120041 ) whose 
parameters are listed in Table |2l 
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Fig. 8. — HD 104237-A Chandra-RETGS spectrum (heg -1,+1, and meg -1,+1 summed; 
binsize=0.02A) integrated over the 145 ks observation. The most prominent hues are labeled. 
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Table 2. Characteristics of HD 104237-A and its spectroscopic companion. 



Parameter HD 104237-A companion Ref. 



d 

M^/Mq 

\og{L^/LQ) 

log(Teff) 

Age 
V sin i 
i 

Porh 



116 pc 

2.25 

2.5 ±0.2 

3.87 

2-5 Myr 

12 ±2 kms-i 

18°tn, 23°t8 

100 ± 5 hrs 

19.859 days 

0.66486 



1.75 

0.5 
3.675 



1 

2 

3 

2 

2 

2,5 

4 

5,3 

3 

2 

2 



Referenc es: (1) iPerrvman et al. 



Bohm et al 



Donati et al. 



(2004); (3) 



Bohm et al 



19971 ): (5) iGradv et all (|2004l ) 



( 19971) : 
(l2006h : 



''Periodicity of modulation observed in Ha, 
possibly associate d to rotational modulation 



Bohm et al. 



(200^ 



(2) 
(4) 

and 

by 



'^Eccentricity of the binary orbit. 



Table 3. Parameters of the HETG observations and count rates of X-ray sources 

HD 104237-A, B, C, D, E. 



Obs ID Start date and time t^xp [ks] 



count rate [cts/ks] 



A 



B 



C 









nth a 


1 st b 


7319 


2006-04-11 23:46:08 


43.39 


23.0 


33.7 


7320 


2006-04-13 10:13:06 


43.83 


19.8 


28.0 


6444 


2006-04-16 09:51:47 


12.57 


20.4 


27.8 


7326 


2006-06-13 05:32:27 


45.61 


35.0 


46.9 



D ^ E 



0.67 0.21 1.5 3.7 

0.62 0.73 1.4 2.3 

1.85 0.48 1.9 3.9 

0.85 0.31 2.2 5.1 



''count rate from zero order counts in the 0.5-8 keV energy range. 

''count rate from first order (MEG -1,+1, HEG -1,-fl) dispersed counts, in the 
1.5-26A wavelength range. 
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These fluxes were determined by fitting each narrow (~ 0.5 A) region of interest with a 
sum of Gaussians and a plasma model continuum, all multiplied by a function of absorption 
by neutral gas with solar abundances. The model was folded through the instrument response 
for each of the negative and positive first orders for each grating, then the positive and 
negative orders were combined before computing the statistic, heg and meg orders were 
combined when feasible (roughly between 3-16 A, where the spectral coverage overlaps with 
significant sensitivity). The line widths were usually frozen at a negligible intrinsic value (in 
§4.3.31 we describe an independent analysis of line profiles). Only for O vill 19 A did we 
need to free the width parameter; this is because of the high dispersion at this wavelength 
and the fact that the H-Lya-like lines are doublets, with a separation of 0.006 A for O Vlll. 

The continuum model was determined by a prior fit of a three temperature-component 
continuum model to relatively line-free regions of the spect rum, using the A strophysical 
Plasma Emission Database continuum emissivities (APEEU; ISmith et al.ll200ll ). This con- 



tinuum then became an a priori determined component for the line fits. For the absorption, 
we used the wabs function from the XSPEC module in ISIS and a fixed value of A'h (see §4.3p . 
From these fits we obtained a single (unabsorbed) flux value as well as a wavelength. 

We did not physically merge any data or response files into new files. Combination of the 
_[_2st Qr(jgj.g -\ya^g clone dynamically. There are many advantages to this technique. Foremost 
is that the individual responses are still applied to the model to compute the predicted 
counts for each order. Second, individual orders can be noticed and ignored as required 
for each region. Third, the combination after folding, and before the statistic, increases 
the signal-to-noise per bin (which does not occur for a joint fit). Fourth, the method is 
transparent to combination of orders from one observation or from multiple observations, 
since all responses, counts spectra, and exposure times are still unique. The visualization in 
Figure [S] shows the summed, combined, and rebinned heg plus meg counts spectrum, but 
this was not the object directly fit. 

Since we have a fairly low-counts spectrum, we also regridded the data (and corre- 
sponding response) to the heg binsize (0.0025 A). This allowed us to combine heg and meg 
spectra, but to keep heg resolution in the heg counts spectrum where feasible. Otherwise, 
we grouped data by an integral numbers of bins independently for each feature. 

For Fexxv (1.85 A), we also included the zeroth order spectrum in the fit, since the 
zeroth order has a comparable number of counts, and though lower resolution, there are no 
significant confusing lines. 



'^APED is available from http://cxc.harvard.edu/atomdb/ 
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These dynamic functions for combining data, regridding responses and data, and for 
grouping data are standard features of ISIS. 
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Table 4. Parameters and uncertainties of 2-T models fitting the ACIS-I spectra of 

HD 104237-A. 



Parameter ^ 


ObsID 2404 


ObsID 3428 


Ti (10^ K) 


4.7 [4.1-5.3] 


4.8 [4.1-5.3] 


EMi (10^2 cni-3) 


19.6 [17.7-20.9] 


22.9 [13.7-26.6] 


T2 (10^ K) 


36.8 [26.8-56.7] 


18.8 [16.3-21.4] 


EMa (10^2 pj^-3) 


3.4 [2.9-4.2] 


15.4 [9.9-17.4] 


L, (1030 erg s^i) 


2.7 


4.1 


Lt (1030 erg s"!) 


3.0 


4.7 


HETGS modeled zero order count rate 


0.020 


0.035 


ACIS-I modeled (unpiled) count rate 


0.16 


0.27 


ACIS-I observed count rate 


0.12 


0.15 



^The ItT confidence intervals are listed in square brackets. Nji is fixed 
at 2 X 10^^ cm~2 and element abundances are fixed to the values derived 
from the hetgs spectrum (see Table [8]). Lg and Lt are the X-ray luminosity 
values in the 0.5-2.0 keV and 0.5-8.0 keV respectively, as in FLG03; these 
value are corrected for absorption. 
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Table 5. Parameters and uncertainties of models fitting the zero-order spectrum of 

HD 104237-A. 



Parameter ^ 


I^ 


II '^ 




2T 


2T 


Nui (1022 em-2) 


0.17 [0.09-0.25] 


0.12 [0.07-0.19] 


Ti (10^ K) 


4.7 [4.3-5.4] 


5.0 [4.3-5.3] 


EMi (10^2 cm-3) 


16.4 [11.6-23.8] 


11.4 [10.6-12.1] 


A^H2 (1022 em-2) 


=iVHl 


=Nhi 


T2 (10*^ K) 


27 [24-31] 


27 [24-30] 


EM2 (10^2 cm-3) 


7.1 [6.0-8.5] 


6.8 [6.3-7.4] 


T3 (10^ K) 


- 


- 


EM3 (10^3 ^j^-3) 


- 


- 


/x (10-12 erg em-2 g-i) 


1.7(1.1) 


1.5(1.1) 


Lx (lO^o erg s-^) 


2.9 (1.7) 


2.4 (1.7) 


0/Oq 


=0.65 


0.46 [0.26-0.63] 


Ne/Neo 


=1.1 


1.7 [1.4-1.9] 


Mg/Mg0 


=0.66 


1.1 [0.6-1.3] 


Si/Si© 


=0.91 


1.2 [0.7-1.5] 


S/So 


=1.24 


1.4 [0.6-1.6] 


Fe/Feo 


=0.50 


0.42 [0.41-0.53] 


Xl (d.o.f.) 


1.12 (149) 


1.05 (143) 



^The 90% eonfidenee intervals are listed in square braekets. 
For X-ray flux and luminosity (/x, Lx; derived in the 0.5- 
7.0 keV energy range) we list the value eorreeted for absorp- 
tion, while the observed (absorb ed) values are in parenth eses. 



Abundances are relative to solar (jAnders fc Grevesselll989l ). 

''Model with two temperature components and abundances 
fixed to the values found from the analysis of the high resolution 
spectra (see i j4.3.1|) . This is our favored model. 

■^Model with two temperature components and variable abun- 
dances. 
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Table 6. Parameters and relative uncertainties for isothermal models fitting the zero-order 

spectra of X-ray sources B, C, D, E. 



Parameter 




B 


C 


D 


E 


Nu (1022 cni-2) 




0.0 [0.0-0.8] 


0.8 [0.-1.7] 


0.44 [0.26-0.67] 


3.4 [2.6-5.2] 


T (10^ K) 




8.4 [1.8-10.] 


7.0 [2.1-32.] 


2.8 [2.0-3.6] 


22. [18-34] 


EM (10^3 pj^-3) 




0.085 [0.064-0.11] 


0.13 [0.012-27.] 


2.9 [0.5-18.] 


1.8 [0.5-3.8] 


/x (10-13 erg cni-2 s" 


-') 


0.35 (0.35) 


0.4 (0.09) 


6.5 (0.73) 


7.1 (1.9) 


Lx (10^0 erg s-^) 




0.057 (0.057) 


0.07 (0.0015) 


1.0 (0.12) 


1.1 (0.31) 


xl 




1.45 


1.24 


1.54 


0.86 


d.o.f. 




13 


2 


6 


20 



Note. — The 90% confidence intervals are listed in square brackets. For X-ray flux and lumi- 
nosity (/x, Lx; derived in the 0.5-7.0 keV energy range) we list the value corrected for absorption, 
while the observed (absorbed) values are in parentheses. For all sources the global metal abun- 
dances arc fixed to 0.3 times solar. 
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Table 7. Line Measurements 





Use'' 




At" 


Ao'' 


f 


C 


h' 


5/g 




Ion 


logT"^ 


Sx'' 


Fexxv 


E 


7.8 


1.861 


1.868 (4.8 


1.15 


0.34 


0.712 


0.441 


1.30 


Ar XVIII 


E 


7.7 


3.734 


3.738 (15. 


0.18 


0.23 


0.154 


0.025 


0.11 


Arxvii 


E 


7.4 


3.949 


3.949 (8.5 


0.56 


0.47 


0.223 


0.335 


0.71 


Arxvii 




7.3 


3.968 


3.953 (15. 


0.001 


0.37 


0.061 


-0.060 


-0.16 


Sxvi 




7.5 


3.992 


4.007 (4.7 


0.26 


0.27 


0.249 


0.006 


0.02 


Sxv 


E 


7.3 


4.088 


4.088 (15. 


0.07 


0.18 


0.064 


0.006 


0.03 


Sxvi 


E 


7.6 


4.730 


4.729 (11. 


1.65 


0.85 


1.221 


0.431 


0.51 


Si XIV 


E 


7.4 


4.947 


4.947 (- 


0.11 


0.19 


0.086 


0.019 


0.10 


Sxv 


E 


7.2 


5.039 


5.040 (2.0 


2.22 


0.61 


1.670 


0.551 


0.91 


Sxv 


E 


7.2 


5.065 


5.065 (- 


0.23 


0.32 


0.390 


-0.160 


-0.49 


Sxv 


E 


7.2 


5.102 


5.107 (5.8 


1.03 


0.46 


0.576 


0.453 


0.99 


Si XIV 


E 


7.4 


5.217 


5.219 (2.8 


0.86 


0.46 


0.407 


0.450 


0.99 


Si XIII 




7.1 


5.285 


5.285 (- 


0.002 


0.18 


0.062 


-0.059 


-0.33 


Si XIII 


E 


7.1 


5.405 


5.405 {- 


0.14 


0.29 


0.136 


0.000 


0.00 


Si XIV 


E 


7.4 


6.183 


6.181 (2.1 


3.51 


0.88 


3.053 


0.462 


0.53 


Si XIV 




7.1 


6.265 


6.256 (15. 


0.45 


0.57 


0.072 


0.374 


0.65 


Mgxii 




7.2 


6.580 


6.580 (- 


0.001 


0.11 


0.063 


-0.062 


-0.59 


Si XIII 


E 


7.0 


6.648 


6.650 (1.1 


3.63 


0.48 


3.609 


0.018 


0.04 


Si XIII 




7.0 


6.687 


6.687 (- 


0.18 


0.21 


0.688 


-0.506 


-2.44 


Si XIII 




6.8 


6.720 


6.720 (- 


0.60 


0.28 


0.035 


0.561 


2.02 


Si XIII 


E 


7.0 


6.740 


6.743 (1.3 


2.06 


0.35 


1.465 


0.594 


1.69 


Mgxii 


E 


7.2 


7.106 


7.107 (5.5 


0.83 


0.32 


0.458 


0.370 


1.15 


Al XIII 




7.4 


7.171 


7.167 (11. 


0.40 


0.28 


0.686 


-0.284 


-1.02 


Mgxi 




6.9 


7.310 


7.317 (6.8 


0.56 


0.31 


0.051 


0.508 


1.65 


Fe XXII 


E 


7.1 


7.681 


7.667 (15. 


0.06 


0.19 


0.037 


0.020 


0.11 


Alxii 


E 


7.0 


7.757 


7.761 (15. 


0.62 


0.37 


0.243 


0.374 


1.00 


Alxii 




6.9 


7.805 


7.790 (15. 


0.06 


0.19 


0.079 


-0.015 


-0.08 


Mgxi 


E 


6.9 


7.850 


7.852 (5.8 


0.57 


0.33 


0.349 


0.218 


0.66 


Alxii 


E 


6.9 


7.872 


7.879 (15. 


0.26 


0.30 


0.220 


0.038 


0.12 


Fexxiii 




7.2 


7.901 


7.897 (15. 


0.33 


0.35 


0.062 


0.267 


0.76 


Fe XXIV 


E 


7.4 


7.986 


7.988 (10. 


0.35 


0.38 


0.192 


0.161 


0.43 


Fe XXIV 


E 


7.4 


7.996 


8.011 (7.5 


0.28 


0.29 


0.097 


0.179 


0.63 


Fe xxiii 


E 


7.2 


8.304 


8.300 (11. 


0.28 


0.33 


0.229 


0.052 


0.16 


Fe XXIV 


E 


7.4 


8.316 


8.317 (15. 


0.20 


0.31 


0.217 


-0.018 


-0.06 


Fe XXIV 




7.4 


8.376 


8.389 (- 


0.34 


0.34 


0.083 


0.253 


0.75 


Mgxii 


E 


7.2 


8.422 


8.422 (1.6 


3.65 


0.50 


3.626 


0.026 


0.05 


Fexxi 


E 


7.1 


8.574 


8.567 (4.9 


0.58 


0.21 


0.188 


0.395 


1.88 


Fexxiii 


E 


7.2 


8.815 


8.816 (3.9 


0.69 


0.28 


0.247 


0.447 


1.59 


Fexxii 


E 


7.1 


8.975 


8.976 (6.6 


0.42 


0.28 


0.321 


0.095 


0.34 


Mgxi 


E 


6.8 


9.169 


9.169 (1.3 


4.24 


0.78 


3.046 


1.195 


1.54 


Fexxi 




7.1 


9.194 


9.194 (- 


1.14 


0.44 


0.144 


0.994 


2.24 


Mgxi 




6.8 


9.230 


9.230 (- 


1.39 


0.50 


0.489 


0.895 


1.80 


Mgxi 




6.8 


9.314 


9.315 (2.7 


2.20 


0.65 


1.447 


0.755 


1.17 


Fe XXII 




7.1 


9.393 


9.393 (- 


0.001 


0.20 


0.072 


-0.071 


-0.36 


Nex 




7.0 


9.481 


9.472 (3.2 


0.48 


0.25 


0.379 


0.096 


0.38 


Fexix 




6.9 


9.695 


9.680 (3.8 


0.72 


0.36 


0.122 


0.601 


1.66 


Nex 


E 


7.0 


9.708 


9.708 (- 


1.35 


0.61 


0.858 


0.492 


0.81 


Fexx 


E 


7.0 


9.727 


9.727 (- 


0.51 


0.39 


0.143 


0.362 


0.94 


Ni XIX 


E 


6.8 


10.110 


10.111 (9.6 


0.31 


0.36 


0.233 


0.074 


0.20 



-27- 



Table 7 — Continued 



Ion 


Use^ 


logrb 


At" 


Ao^ 


Fexx 


E 


7.0 


10.120 


10.127 (14.) 


Nex 


E 


7.0 


10.239 


10.240 (1.7) 


Fe XXIV 


E 


7.4 


10.619 


10.619 (-) 


Fexix 


E 


6.9 


10.641 


10.641 (-) 


Fexix 


E 


6.9 


10.649 


10.649 {-) 


Fe XXIV 


E 


7.4 


10.663 


10.663 (-) 


Fe XXIII 




7.2 


10.981 


10.983 (13.) 


Neix 


E 


6.7 


11.001 


11.001 (-) 


Fe XXIII 


E 


7.2 


11.019 


11.019 (-) 


Fe XXIV 


E 


7.4 


11.029 


11.029 (-) 


Fe XVII 




6.7 


11.131 


11.123 (15.) 


Fe XXIV 


E 


7.4 


11.176 


11.176 (2.3) 


Fe XVII 


E 


6.7 


11.254 


11.251 (2.8) 


Fe XVIII 


E 


6.8 


11.326 


11.321 (4.5) 


Fe XVIII 


E 


6.8 


11.527 


11.524 (6.8) 


Neix 


E 


6.6 


11.544 


11.548 (4.9) 


Fe XXIII 




7.2 


11.736 


11.741 (1.8) 


Fe XXII 


E 


7.1 


11.770 


11.773 (2.3) 


Nex 


E 


6.9 


12.135 


12.131 (1.0) 


Fe XXIII 




7.2 


12.161 


12.150 (4.2) 


Fe XVII 


E 


6.7 


12.266 


12.265 (3.4) 


Fexxi 


E 


7.1 


12.284 


12.284 (2.1) 


Fexx 


E 


7.0 


13.385 


13.396 (6.5) 


Fexix 


E 


6.9 


13.423 


13.423 (-) 


Neix 


E 


6.6 


13.447 


13.449 (1.5) 


Fexix 


E 


6.9 


13.462 


13.462 (-) 


Fexix 




6.9 


13.497 


13.497 {-) 


Fexix 


E 


6.9 


13.518 


13.519 (1.9) 


Neix 




6.6 


13.552 


13.554 (3.3) 


Fexix 


E 


6.9 


13.645 


13.653 (-) 


Neix 


E 


6.6 


13.699 


13.694 (6.4) 


Fexx 


E 


7.0 


13.767 


13.762 (10.) 


Fexix 


E 


6.9 


13.795 


13.790 (7.6) 


Fe XVII 


E 


6.7 


13.825 


13.827 (11.) 


Fe XVIII 


E 


6.8 


14.208 


14.207 (8.3) 


Fe XVIII 


E 


6.8 


14.256 


14.249 (13.) 


Fexx 


E 


7.0 


14.267 


14.268 (13.) 


Fe XVIII 


E 


6.8 


14.343 


14.355 (15.) 


Fe XVIII 


E 


6.8 


14.373 


14.374 (5.4) 


Fe XVIII 




6.8 


14.425 


14.410 (15.) 


Fe XVIII 




6.8 


14.534 


14.547 (6.1) 


Fexix 


E 


6.9 


14.664 


14.667 (11.) 


Oviii 


E 


6.7 


14.821 


14.817 (4.4) 


Fe XVII 


E 


6.7 


15.014 


15.014 (3.4) 


Fexix 


E 


6.9 


15.079 


15.087 (9.2) 


Oviii 


E 


6.7 


15.176 


15.178 (15.) 


Fexix 


E 


6.9 


15.198 


15.209 (11.) 


Fe XVII 


E 


6.7 


15.261 


15.261 (5.5) 


Fe XVII 


E 


6.7 


15.453 


15.451 (15.) 



ft' 



5p 



Sx" 



0.39 


0.27) 


3.98 


0.83) 


1.36 


0.60) 


1.05 


0.78) 


0.42 


0.60) 


0.80 


0.58) 


0.39 


0.47) 


1.23 


0.61) 


1.34 


0.72) 


1.25 


0.65) 


0.001 


0.09) 


2.72 


0.77) 


2.54 


0.76) 


1.61 


0.64) 


0.99 


0.61) 


1.67 


0.69) 


4.96 


0.99) 


3.34 


1.01) 


24.0 


(2.9) 


4.88 


2.05) 


3.60 


1.37) 


7.82 


1.75) 


1.53 


1.17) 


2.57 


1.38) 


15.2 


(2.7) 


2.34 


1.73) 


0.36 


0.74) 


9.60 


2.34) 


3.95 


1.57) 


1.49 


2.46) 


6.91 


3.20) 


3.77 


3.37) 


3.90 


3.33) 


3.92 


3.26) 


9.60 


4.83) 


3.57 


3.68) 


3.59 


4.22) 


1.52 


3.24) 


9.40 


3.77) 


3.49 


2.96) 


4.36 


3.49) 


4.95 


3.14) 


4.37 


1.78) 


32.8 


(7.4) 


6.55 


6.59) 


3.26 


6.12) 


7.39 


6.26) 


17.3 


(6.8) 


4.26 


4.53) 



0.161 


0.232 


0.87 


2.817 


1.165 


1.41 


1.576 


-0.217 


-0.36 


0.273 


0.772 


1.00 


0.264 


0.152 


0.25 


0.830 


-0.026 


-0.04 


1.443 


-1.053 


-2.26 


0.418 


0.812 


1.33 


0.953 


0.382 


0.53 


1.047 


0.205 


0.32 


0.810 


-0.809 


-8.82 


1.915 


0.799 


1.04 


1.146 


1.397 


1.84 


0.838 


0.766 


1.19 


0.890 


0.101 


0.17 


1.362 


0.304 


0.44 


3.323 


1.637 


1.66 


3.561 


-0.223 


-0.22 


24.246 


-0.254 


-0.09 


1.934 


2.943 


1.43 


3.113 


0.483 


0.35 


7.538 


0.282 


0.16 


1.043 


0.487 


0.42 


0.559 


2.008 


1.46 


12.577 


2.657 


1.00 


1.274 


1.066 


0.62 


2.245 


-1.884 


-2.54 


4.958 


4.646 


1.98 


1.904 


2.041 


1.30 


0.801 


0.687 


0.28 


6.299 


0.607 


0.19 


0.917 


2.856 


0.85 


2.036 


1.863 


0.56 


2.787 


1.129 


0.35 


12.652 


-3.051 


-0.63 


2.445 


1.122 


0.31 


1.561 


2.030 


0.48 


1.429 


0.090 


0.03 


3.077 


6.324 


1.68 


0.668 


2.821 


0.95 


2.365 


1.999 


0.57 


1.558 


3.389 


1.08 


1.174 


3.193 


1.79 


39.169 


-6.377 


-0.86 


1.809 


4.738 


0.72 


2.741 


0.515 


0.08 


1.514 


5.875 


0.94 


11.351 


5.939 


0.87 


1.484 


2.778 


0.61 
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4-3.1. Emission Measure Distribution and Abundances 

We reconstructed an emission measure distribution (EMD) and elemental abundances 
by minimizing line flux residuals for features marked in Table [7] with an "E" in the "Use" 
column. This was an iterative process, since some features are blended or mis-identified: we 
excluded those with very large residuals during the fitting process. We also excluded the 
He-like resonance and intercombination lines since they are density dependent. Theoreti- 
cal emissivities from APED were integrated over the temperature range of sensitivity after 
weighting by the trial EMD and abundances. The fit was subject to the constraint that the 
EMD be smooth. To obtain some estimate of the uncertainty in the reconstructed EMD 
and abundances, we performed 100 Monte-Carlo iterations in which the observed fluxes were 
randomly perturbed using a Gaussian distribution with their measurement uncertainties as 
the Gaussian a. We computed the standard deviation in the resulting EMD at each tem- 
perature and for the abundances and use these for the fit uncertainty. This is not rigorously 
correct because it does not account for correlations in the allowed solutions. Other ex amples 



of this method, and further details, can be found in iHuenemoerder et al.l (120061 . 120071 ). along 
with requisite caveats and limitations of emission measure modeling. As a final check, the 
synthetic spectrum was evaluated and compared to the observation (see Figure [TOl) . Model 
line fluxes are listed in Table [71 the EMD is shown in Figures [9], and abundances are given 
in Table El 

The emission measure distribution shown in Figure [9] is characterized by significant 
amount of plasma over a wide temperature range, as anticipated on the basis of a visual 
inspection of the spectrum, and by a peak at logT[K] ~ 6.9. These characteristics are 
simi lar to the temperature di stributions derived for several coronae of late-type stars (see 
e.g.. lSanz-Forcada et al.ll2003l ). typically presenting an emission measure peak around 10^ K. 
However, coronae of active stars with X-ray luminosity > lO'^^erg s~^ often present an addi- 
tional significant peak at higher temperatures (a few 10^ K) typically associated with flaring 
activity. Hot plasma is present in HD 104237-A as well, but with a lower relative weight 
than observed in other stellar coronae at a similar activity level. 

In Figure [9] we also compare the EMD model with the 2T model obtained from the 
low-resolution spectra (model I in Table [5]): the emission measure values for the components 
of the 2T model are scaled for comparison purposes, but their sum compares well with 
the integrated emission measure distribution. Figure [101 shows the synthetic spectrum from 
the EMD model, compared with the actual spectrum. Figure [11] compares the predicted 
spectrum from the two models with the actual observed spectrum and shows that while 
the overall characteristics of the 2T model are somewhat similar to the detailed EMD 
model, the latter reproduces the spectral features of the high-resolution spectrum much 
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Table 7 — Continued 



Ion 



Use'' logT^ \f^ 



Xo^ 



ft' 



sp 



Sx"" 



Fe XVIII 




6.8 


15.494 


15.488 


(15. 


2.17 


3.32) 


0.386 


1.784 


0.54 


Fe XVIII 


E 


6.8 


15.625 


15.625 


(15. 


2.51 


3.69) 


3.604 


-1.096 


-0.30 


Fe XVIII 




6.8 


15.759 


15.763 


(15. 


3.37 


4.31) 


0.497 


2.876 


0.67 


Fe XVIII 


E 


6.8 


15.824 


15.826 


(15. 


4.60 


5.04) 


2.251 


2.346 


0.47 


Fe XVIII 


E 


6.8 


15.870 


15.875 


(15. 


5.22 


5.49) 


1.200 


4.018 


0.73 


Oviii 




6.7 


16.006 


16.004 


(2.2 


21.7 


(5.0) 


9.436 


12.242 


2.43 


Fe XVIII 


E 


6.8 


16.071 


16.075 


(6.1 


4.22 


2.39) 


5.304 


-1.087 


-0.45 


Fexix 




6.9 


16.110 


16.095 


(2.5 


7.60 


2.97) 


2.642 


4.963 


1.67 


Fe XVIII 




6.8 


16.159 


16.150 


(4.0 


6.05 


2.49) 


2.132 


3.922 


1.57 


Fe XVII 


E 


6.7 


16.780 


16.779 


(4.0 


41.4 


11.6) 


21.558 


19.810 


1.71 


Fe XVII 


E 


6.7 


17.051 


17.052 


(6.3 


29.4 


10.9) 


26.525 


2.838 


0.26 


Fe XVII 


E 


6.7 


17.096 


17.100 


(4.4 


37.8 


11.2) 


24.814 


13.014 


1.16 


Fe XVIII 




6.8 


17.623 


17.635 


(4.7 


6.14 


3.39) 


4.127 


2.010 


0.59 


Ovii 




6.4 


17.768 


17.762 


(8.1 


4.87 


3.26) 


0.663 


4.203 


1.29 


Fe XVII 


E 


6.6 


17.793 


17.778 


(15. 


0.001 


1.48) 


0.000 


0.000 


0.00 


Ovii 




6.4 


18.627 


18.627 (- 


0.001 


1.24) 


2.329 


-2.328 


-1.87 


Ca XVIII 


E 


7.1 


18.691 


18.691 (- 


0.46 


3.05) 


0.644 


-0.185 


-0.06 


O viii 


E 


6.7 


18.970 


18.971 


(2.2 


104. 


(17.) 


89.094 


14.826 


0.85 


Caxvi 




6.7 


20.859 


20.869 


(12. 


7.09 


7.70) 


0.043 


7.046 


0.92 


Nvii 




6.5 


20.910 


20.881 


(30. 


0.001 


5.28) 


1.958 


-1.958 


-0.37 


Ovii 


E 


6.3 


21.601 


21.605 


(10. 


20.0 


16.3) 


28.511 


-8.519 


-0.52 


Ovii 




6.3 


21.802 


21.817 


(6.4 


17.3 


15.3) 


3.861 


13.392 


0.88 


Ovii 




6.3 


22.098 


22.083 


(15. 


2.20 


9.83) 


16.800 


-14.603 


-1.49 


Nvii 


E 


6.5 


24.782 


24.797 


(15. 


20.4 


14.3) 


13.131 


7.253 


0.51 



''The label E indicates the hnes used for the emission measure reconstruction. 

''Average logarithmic temperature (kelvins) of formation. 

'^Theoretical wavelengths of identification (from APED), in A units. If the line is a multiplet, 
wc list the wavelength of the strongest component. 

'^Measured wavelength, in angstrom units (uncertainty is in mA). Errors equal to (-) indicate 
that the line wavelength was frozen because the feature is very weak. 

"Measured line fluxes in 10^® photons cm^^ s^' for HEG-I-MEG first orders, with la uncer- 
tainties, assuming Nii= 1.5 X 10^^ cm~'^. 

Fluxes (in 10~® photons cm~^ s~^) predicted using the emission measure distribution model 
(see Fig. [gjl. 

^Linc flux residuals, <5/ = /; — /t. 

''Discrepancy between observed and predicted fiuxes in units of a. 



more satisfactorily. 
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Fig. 9. — Emission measure distribution inferred from the measured line fluxes. Crosses 
show the EM component derived from the 2T fit to zero-order spectra. Since we are plotting 
the emission measure integrated over the temperature bins, we applied an arbitrary scaling 
factor of 0.1 to the 2T models in order to compare the relative weights with the shape of the 
EMD; total volume emission measures for the two models, 2T and EMD, are comparable, 
as they should be. 



The element abundances of the X-ray emitting plasma of HD 104237-A derived from 
the analysis of its hetg s spectrum are listed in Table [HI both compared to solar values 
( Greyesse fc Sauvallll998l ) and to the photospheric values of HD 104237-A from lAcke &: Waelkens 



(120041 ). These abundances show a general depletion of elements in the X-ray emitting plasma 
with the exception of S, generally in line with X-ray studies of other PMS stars (e.g.. 



Maggio et al.l 120071 ). however with some possible differences: the Ne/Fe abundance ratio 
is relatively lo w (^ 2), with respect to the extremely hig h values observed in other young 



stars (see e.g., iKastner et al.l l2002l : lArgiroffi et al.l l2005l finding values of Ne/Fe ~ 10 for 



PMS stars in TWA); also, possible dependence of fractionation on the element first ion iza- 
tion potential, as often observed in coronae of late- type stars (e.g.. iTelleschi et al.ll2005l ). is 
not obviously present in the X-ray emitting plasma of HD 104237-A. 
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Fig. 10. — Comparison of the heg+meg spectrum of HD 104237-A with the model (hght 
histogram) assuming A^h= 1-5 x 10^^ cm~^. Residuals are also shown. 
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Table 8. Element abundances of X-ray emitting plasma from high resolution spectrum. 



Element 



AIA^' 



^Mphot*^ 



^photMo [dex] 






0.65 ±0.14 


0.82 ±0.33 


-0.10 ±0.14 


Ne 


1.05 ±0.10 






Mg 


0.66 ±0.08 


1.02 ±0.20 


-0.19 ±0.07 


Si 


0.91 ±0.10 


0.87 ±0.25/0.45 ±0.11 


0.02 ±0.11/0.31 ±0.09 


S 


1.24 ±0.36 


1.75 ±0.80 


-0.15 ±0.14 


Fe 


0.50 ±0.05 


0.42 ±0.19 


0.08 ±0.18 



^Element abundances of the X-ray emitting p lasma from the high resol ution 
X-ray spectrum, expressed as relative to solar (JGrevesse fc Sauvallll998l ). We 
also list Ict errors. 

''Element abundances of the X-ray emitting plasma from the high resolution 
X-ray spectrum, expressed as relative to photospheric abundances. Errors are 
calculated propagating the errors on both the photospheric abundances, and 
the abundances of the X-ray emitting plasma. 



■^Photospheric values of clement abundances from lAcke fc Waelkend (J2004( ): 
for Si we list their values derived from Sii, and from Sin respectively. 




Fig. 11. — Comparison of the heg+meg spectrum of HD 104237-A with the spectrum 
simulated using the EMD model {left) and the 2T model [right). Note the different scales 
onx- 
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4-3.2. Density and Temperature Diagnostics from Line Ratios 



The He-like triplets, in particular the cooler Ne ix and O Vll triplets, t hanks to their 



sensitivity to plasma density and temperature (e.g., iGabriel fc Jordan! Il969l ). and to pho- 



toexcitation from possible UV f ield, provide effectiv e diagnostics for the X-ray emission 
mechanism in young stars (e.g., iKastner et al.ll2002l ). In particular, in absence of intense 
UV field, the ratio of the forbidden to intercombination line, R = f/i, is sensitive to plasma 
densities above a critical density value that depends on the element and increases from ap- 
proximately 3 X 10^° cm~^ for Ovil to about 3 x 10^^ cm~^ for Sixill. The presence of 
strong UV field, either photospheric for hot stars or produced in accretion shocks in young 
stars, affects the R ratio similarly to high densities, i.e. decreasing its value. In the case 
of HD 10423 7- A the stellar effective temperature is low enough (see Table |2]) to rule out 
any effect of the photospheric UV field on the X-ray triplet lines. The ratio of the summed 
forbidden and intercombination lines with the resonance line, G = {f + i)/r, provides instead 
a temperature di agnostic (though w ith some caveats on the reliability of the diagnostics, as 
discussed e.g., by lTesta et al.ll2004l ). 



The O VII He-like triplet at 22A is too weak in hetgs spectrum of HD 104237-A to 
provide any meaningful diagnostic (only 3 counts are detected for the strongest triplet line, 
the 21. 6A resonance line). Neix, Mgxi, and Sixill have good signal-to-noise ratios and 
provide useful info on the plasma density and temperature. The Sixill triplet is at the 
low density limit as for all other observed high resolution spectra of late-type stars (e.g.. 
Testa et al.l[20oi ). 



In Figure [12] we show the Ne ix and Mg xi triplet regions with our best fit parametric 
models (upper panels), and the confidence contours in the parameter space of the line ratios R 
and G (lower panels). In the R and G p arameter spa c e we overplot the grid of corresponding 
densities and temperature according to lSmith et al.l (1200 ll ). The fitting model for the Neix 
triplet re gion includes the blends o f mostly Fe lines as computed from the EMD model 



Huenemoerder et al.ll2006l ). The best fit to the Neix triplet lines corresponds to a 



[see e.g., 

density rie ~ 6 x 10^^ cm~^; however the relatively low statistics of the spectrum does not 
allow to put stringent limits on this value and we can reliably derive only an upper limit to 
the density of about ne < 6 x 10^ 



U2 pj^ 3_ 



The fits to the Mgxi t riplet region included Ne x Ly series blends, and blends of Fe lines 



(see e.g.. 



Testa et al.ll2004l ). Mgxi allows to put slightly better constraints on the plasma 
density: the best fit density is ne ~ 7 x 10^^ cm^'^ with a la interval 1 — 18 x 10^^ cm~^, 



somewhat high compared with typical values of active late-type stars (e.g.. lTesta et al.ll2004l ). 
Also the G ratio is somewhat extreme when co mpared with typic al coronal values, being at 
the low-end of the range of measured ratios (JTesta et al.l |2004| ). even though the rather 
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large uncertainties do not allow to draw definite conclusions. Interpreting the G ratio as a 
temperature diagnostics would point to a plasma temperature in agreement with the expected 
values for the formation of these lines. This is at variance with findings relative to He-like 
triplets of stellar coronae whose G ratios generally underestimate the plasma temperature. 
Therefore it is difficult to interpret an unusually low G ratio in HD 104237-A as a temperature 
effect considering that its temperature structure does not seem peculiar when compared to 
stellar coronae. We note that also the Neix G ratio has a v alue among the l owest measured 
in a large sample of active stars at different activity levels ( jNess et al.ll2004l ). 



I 



fTI 



mw 




13.55 13.5 

Wavelength (A) 



Wavelength (A) 




Fig. 12. — Chandra-KE.TG'S, spectrum of HD 104237-A in the Neix {left) and in the Mgxi 
[right) He-like triplet wavelength regions (black; rebinned by a factor 2 with respect to 
original wavelength grid) with best fit model (red). The lower panels show the corresponding 
confidence contours (68%. 90%, 98%) in the R — G parameter space and indicate the density 
and temperature diagnostics derived from the triplet line ratios. 
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4-3.3. Line profiles 

Spectral line profiles in X-rays allow us to probe the bulk velocity of hot plasma and 
its turbulent state. We carried out the measurement of line centroid positions of several 
prominent emission lines detected in the meg spectra with sufficient counts for meaningful 

fits. 

The methodology used here is somewhat different from the line-flux measurements de- 
scribed in §4.31 Line centroids are very sensitive to blending, even with weak features. Here 
we take advantage of a model spectrum which predicts the relative positions and relative 
strengths of many features in any region. Including thermal broadening in the flux evaluation 
automatically accounts for any line-width difference between species, and turbulent broad- 
ening can be included to test for excess width. For the best accuracy, these flts were only 
done on strong-line regions, whereas flux measures for emission measure reconstruction were 
also done for weak lines, since the absence of a line can be signiflcant for the temperature 
structure. Since the line widths are barely resolved, if at all, the width measurement can 
best be done on only the highest signal-to-noise features, and though the heg has higher 
spectral resolution, the higher counts in the MEG can makes its determination of width more 
signiflcant. 

Utilizing the latest calibration product for MEG line spread function (available for Chan- 
dra CALDB v.3.4.0 or higher), grating ARFs and RMFs are generated for meg +1 and -1 
orders. Then the ±1 order meg spectra are simultaneously flt with a simple 2T model 
similar to the ones listed in Table [5] (close to model I) with line centroid and broadening 
parameters as free parameters. We have taken thermal broadening into account upon fltting. 
Since the observations were taken at four different times over the two-month period, each 
spectrum was corrected for appropriate barycentric motion of the Earth (-8.3, -9.9, -9.4, and 
+13 km s~^ for ObsID 6444, 7319, 7320 and 7326, respectively) prior to the fltting process. 

We selected eight prominent lines in a wavelength range from 6 to 20A: the H-like Lya 
transitions of Si, Mg, Ne, and O, and the Fexvil transitions at ~ 15.0 and 17. lA. The 
measured line velocities are tabulated in Table [91 

The mean line centroid velocity of the source is 46 ± 42 km s~^ with a 99% confldence 
level. The large uncertainty in the measurement results from the low S/N level detected 
in each emission lines. As noted, HD 104237-A is a binary system. Generally a small 
radial velocity correction should be made - prior to line fltting - for a spectrum obtained 
at each epoch. However, we intentionally neglect this since (i) we do not know the relative 
contribution of each component to the spectra and (ii) all of the spectra were obtained in 
the phase ranging from 0.25 - 0.5 (see bottom panel of Figure H]), placing the stars near 
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apastron where the change in radial velocity happens to be very small, of order of 7 km s 



based on Bohm et al. 2004 



At the observed times, the radial velocity of the primary star, the HAe star, should 
be observed at ~ +20 km s^^, whereas the companion, K3, should be at ~ +10 km s~^. 
While the measured radial velocity prefers the primary star to be the source of hard X- 
rays, the companion is not statistically ruled out as the source. Furthermore, with the 
radial velocity correction of the systemic binary motion, the measured line centroid velocity 
becomes consistent with the hard X-ray source being at rest. 

We also measured the width of the lines, using meg ±1 (and heg ±1 for lines with 
sufficient counts like the Nex Lya transition) order spectra. The line widths are consistent 
with zero broadening for all the strongest lines, with upper limit for Wturb of the order of 
300 km s~^. The fit to the Fexvil 15A line provides a best fit width corresponding to a 
turbulent velocity fturb of about 280 km s~\ and barely compatible with at the 99% level. 
However, we note that the statistics are limited, and a small error in the continuum can 
cause fturb to be compatible with at the 68% confidence level. Also, the other Fexvil lines 
around 17A are compatible with zero broadening casting doubts on the results to the fit to 
the ISA line. 



5. Discussion 

The HD 104237 young association, with presumably coeval stars spanning O.15-2.3M0, 
offers a laboratory for studying X-ray emission in young low/intermediate mass stars. The 
brightest X-ray source (and also so in optical, UV, and IR bands) in the stellar group is a 
system harboring a binary of intermediate mass stars (2.25/1.75 Mq) with a HAe star, and 
a slightly lower mass star in an early evolutionary stage with spectral type K3. The Herbig 
Ae star HD 104237-A with its proximity, high Lx, and low line-of-sight extinction is one of 
the few HAe stars accessible to high resolution X-ray spectroscopy, allowing us to investigate 
the nature of the X-ray emission from young intermediate mass stars that is still a mystery. 



The age of the binary system is estimated between 2 and 5 Myr (JBohm et al.ll2004l . FLG03) 



and the HAe star is correspondingly in the late stages of its pre-main sequence phase. 

The long Chandra-tiETGS observation provides important clues to the X-ray emission 
mechanism through spectral diagnostics and variability. Different scenarios may be consid- 
ered to explain the X-ray emission of HD 104237-A: (1) all X-ray emission comes from the 
HAe star; (2) all X-ray emission comes from the companion; (3) both binary components 
contribute to the observed X-ray emission. The companion hypothesis is still considered a 
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viable explanation for the X-ray emission of Herbig Ae stars, considering that these inter- 
mediate mass stars are not expected to produce X-rays either through wind shocks, as in 
earher type stars, because their winds are weak, or through magnetic activity as in later type 
stars, because they should not have an efficient solar-like dynamo mechanism producing mag- 
netic fields since they lack sizable convective zones. However, some alternative mechanisms 
might be at work in the early stages of their evolution: (i) X-ray emission from accretion 
shocks, as in some classical T Taur i stars fe.g.jKastner et al.ll2002l) . (ii) coronal activity due 



to non-solar dynamo mech anism (JTout &: Pringk 



19951 ). (iii) magnetically confined winds 



( iBabel fc Montmerldll997l ). All these mechanisms require the presence of some sizable mag- 
netic field in this evolutionary stage. The magnetic field in HAe stars has been investigated 
on large sample of these intermediate mass stars, and a small number of them do show evi- 
dence of r aagnetic fields of strength from a few tens of Gauss to a few hundreds of Gauss in a 
few cases JPonati et al.lll997l : iHubrig et allbood : IWade et al.lboosl . I2OO7I : iHubrig et allboO?! : 
Catala et al.l 120071 ). HD 104237-A is one of those showing that, whatever the mechanism 
that generates it, at some stage these intermediate mass stars can have magnetic fields. The 
expected X-ray properties are different depending on the mechanisms with (i) characterized 
by soft X-rays emitted by high density plasma, and harder X-ray emission due to magnetic 
confinement (and presumably heating) in (ii), and (iii). In particular in the scenario of 
coronal activity we expect a large degree of variability and fiare-like dynamic events. 

What do we know from observations of other HAe stars? X-ray observations with 
Chandra allow to resolve the Herbig Ae stars from close (> 1") compai iions and therefore 
reduce the confusion, and still indicates that a significant fraction (~ 35%; IStelzer et al.ll2006f) 



of Her big Ae star with unknown companions are X-ray sources. As shown by IStelzer et al. 



(J2006|), the HAe stars have X-ray properties (Lx, T) very similar to lower mass young stars; 
therefore they conclude that either their emission and, likely. X-ray production mechanisms 
are very similar to TTS (i.e., coronal), or their X-ray emission comes from yet to be discovered 
companions. A few exceptions e xist: HD 163296 ha s been observed with Chandra and shows 
an unusually soft spectrum, that lSwartz et al.l (120051 ) tentatively attri bute to accretion shoc ks 
as in some CTTS; analogously AB Aur shows a very soft spectrum (JTelleschi et al.ll2007l ). 



The spectroscopic companion to HD 104237-A is a later spectral type s tar (similar to 



SU A ur) and similar stars a re known to be strong X-ray emitters (see e.g., iGetman et al. 



20051 . and lGiidel et al.l 120071 ). Therefore the K3 companion is expected to contribute signif- 
icantly to the overall X-ray emission. We explore the companion hypothesis assuming that 
the X-ray emission is entirely due to the K-type companion to HD 104237-A, and compar- 
ing its properties with the X-ray properties of stars of similar mas s and evolutionary stage. 
The Chandra Orion Ultradeep Project (COUP; iGetman et al.ll2005l ) provides information for 
putting the X-ray emission properties of the HD 104237 members into context. In Figure 
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we coin pare some spectral pro perties of the H P 104237 inembers wi th the large sample from 
COUP fjPreibisch et al.ll2005l ). In particular iPreibisch et al.l (120051 ) analyzed the Chandra- 
ACIS spectra of the COUP sources with a two-temperature model and studied how the best 
fit temperatures change with spectral type. We superimpose on the results from COUP our 
temperature values from the modeling of the zero-order spectra, assuming a spectral type 
K3 for HD 104237-A (Figure [T3| left panel). We analyze in the same context also the results 
found for source D and E, which have well determined spectral types. 
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Fig. 13. — Comparison of plasma temperature of the HD 104237 members, obtained from 2- 
T (for source A) or 1-T (f or source D, and E) an alysis of their zero-ord er spectra, with results 
from the COUP sample (IPreibisch et al.ll2005l ). We used plots from IPreibisch et al.l (120051 ) 
(their Fig. 12, and 13; reproduced by permission of the AAS) and superimposed the values 
found for HD 104237-A, D and E (star symbols of different colors: A = yellow, D = red, and 
E = blue). Left panel: temperature of X-ray emitting plasma vs. effective temperature. Txi 
{crosses), and Tx2 (diamonds) are the two temperatu re values of a 2-T mod el fit to the low 



resolution X-ray spectra. Right panel: Txi vs. Tx2- IPreibisch et al.1 (120051 ) (see references 
therein) include in the plot results for a sample of G- and K-type MS stars [open squares), 
and typical values for structures in the solar corona (BG = background corona; AR = active 
region; FL = flares). 

The plots of Figure [T^ show that the hotter of the two temperatures fitting the spectrum 
of source A compares well with the high temperature component found for early K-type stars 
in the young star forming region of Orion, and also with the temperature of source E which 
has same spectral type as the spectroscopic binary companion to A. On the other hand, the 
soft component seems significantly cooler than the soft component of these similar stars. 
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This soft component is re miniscent of the so f t spectra found for t he two Herbig Ae stars, 
HD 163296 and AB Aur JSwartz et all boosi : IXelleschi et allboOTh . These recent findings 
might lend support to the idea that characteristic soft X-ray emission could be a defining 
property of X-ray spectra of HAe stars. The persistent hot component, however, cannot 
be explained by either the wind models or magnetic infall scenario proposed for the CTTS, 
and therefore the detection of the hot component suggests the presence of magnetically 
confined plasma. This hot emission is likely produced, at least partially, by the later type 
companion; on the other hand the X-ray emission mechanisms - emission from magnetically 
confined coronal plasmas, or from accretion streams, or from magnetically confined winds - 
may not be mutually exclusive and may all contribute to the X-ray emission in some HAe 
systems. Possible generation of magnetic fields in HD 104237-A cannot be easily explained 
with existing models. The estimated position in the HR diagram and evolutionary tracks 
imply that HD 104237-A is already completely radiative and therefore it should not be 
able to produce any ma gnetic field by s olar-l ike dynamo mechanism. In the alternative 
shear magnetic model by lTout &: Pringld (119951 ). X-ray luminosities at levels high enough to 
explain the observed Lx of HAe stars can be produced for time scales of the order of 1 Myr, 
and then they are rapidly decreasing as Ly(t) = Lxr. i } + t/ tn)~^, where Lxq depends on the 
stellar parameters. As calculated by ISkinner et al.l (l2004l ) depending on the set of stellar 
parameters adopted in the ranges discussed in ^ this model underestimates the observed 
Lx by a factor ~ 4 (for spectral type A4ive and age t ~ 2 Myr) to ~ 25 (for the more 
recently determined spectral type A7.5ve-A8ve and age 2-5 Myr). 

HD 163296, AB Aur and HD 104237- A show similar level of X-r ay emission with frac - 
tional luminosity logLx/i^boi of about -5.5 (ISwartz et al.ll2005l ). -5.6 (JTelleschi et al.l 120071 ). 
and -5 (this work; logLx/i^boi~ —4 if we assume that the K3 companion is the X-ray source) 
respectively. Another similarity between HD 104237-A and AB Aur consists in their X- 
ray variability: both sources are slowly varying within a factor ~ 2, do not present any 
clear evidence of flaring activity, and their X-ray lightcurve is compatible with their respec- 
tive rotational period as derived from periodicity of optical lines. The hetgs observations 
presented in this paper do not provide enough phase coverage to definitely determine the 
presence of X-ray modulation on the timescales of modulation of the Ha line, however, if 
the observed modulation is confirmed to be consistent with the rotation period of the HAe 
star, this would imply that a significant portion of the observed X-ray emission is produced 
by the HAe star. X-ray rotational modulation has been previously observ ed in a number 
of X-ray bright PMS stars, such as for instance by iFlaccomio et al.l (120051 ) who found ev- 
idence of X-ray emission being modulated by rotation in at least 10% of the Orion PMS 
stars sample observed in the COUP campaign; these findings are also compatible with theo- 
retical modeling of TTS coronae from extrapolations of surface magnetograms derived from 
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Zeeman-Doppler imaging (JGregory et al.ll2006l ) 



The high resolution hetgs spectra offer additional spectral diagnostics providing clues 
to the yet unestablished X-ray production mechanisms in young intermediate mass stars. 
Particularly useful diagnostics are provided by the He-like triplets: specifically the Neix 
and Mg xi triplets, since the O vil triplet is essentially undetected, and the Si xill lines 
are at the low density limit (therefore providing only an upper limit to density of a few 
10^^ cm~^). Both Ne and Mg He-like triplets line ratios indicate plasma densities somewhat 
larger than typically observed in stellar coronae at the same level of activity, though the 
limited statistics do not provide very stringent constraints to these findings. In the scenario 
tentatively drawn above, with the Herbig Ae star being the main contributor to the cool 
(< 10^ K) emission and the K3 companion producing most of the hard emission, these 
triplets lines would be produced mainly by the Herbig Ae star. Considering the presence 
of significant ongoing accretion (at an estimated rate of ~ 10~^Mq yr~^), in analogy with 
TW Hya an d other accreting T Tauri stars showing high density in their high resolution 



spectra fe.g..lKastner et al.ll2002l : IStelzer &: Schmittll2004l : ISchmitt et al.ll2005l : iGiinther et al. 



20061 : lArgiroffi et al.l 120071 ). the high densities observed for HD 104237-A can be interpreted 
as a signature of X-ray production in accretion shocks. The high resolution spectra also 
provide line profile diagnostics potentially able to determi ne the contribution to the X-ray 
emission of each component of a binary system (see e.g., Ilshibashi et al.ll2006l ). We find 
that the measured line centroids for the strongest lines (see §4.3.3^ indicate that the Herbig 
Ae star is more likely, with respect to the K3 companion, to be the primary X-ray source. 
However, in our case the involved velocities are small, and the statistics limited, therefore not 
providing conclusive results. As discussed in §4.3.3l no convincing evidence of line broadening 
is found in the strong X-ray lines. 

Among the other sources in the association, source D shows interesting characteristics: 
as clear from Figure [T3] (left panel), source D is characterized by an unusually cool plasma 
temperature with respect to the bulk of similar M-type young stars in Orion. Even though 
no clear evidence of ongoing accretion is available for this star, its Ha emission is near the 
CTTS limit indicating that this star might still be actively accreting from its circumstellar 
disk; in this scenario, the soft excess might indicate that X-rays produced by shocked ac- 
creting plasr na dominate the X-r ay emission of this star, analogously to the peculiar case of 
TW Hydrae (JKastner et al.ll2002l ). and differently from the majority of TTS where the cor o- 
nal component generally dominates the overall X-ray emission (e.g.. iPreibisch et al.ll2005l ). 
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Conclusions 



We have analyzed Chandra-UETG observations of young low-mass and intermediate- 
mass stars in the small stellar group associated with the Herbig Ae star HD 104237 to 
investigate the X-ray emission mechanisms in their pre-main sequence phase. The close 
(~ 15 AU) binary system of the HAe star and its K-type companion cannot be resolved in 
X-rays, however our findings with respect to the X-ray variability and spectral properties 
provide clues to the origin of the X-ray emission and the X-ray production mechanisms. The 
modulation on time scales of the estimated rotation period of the HAe star suggests that 
the primary contributes significantly to the observed X-ray emission. The spectral analysis 
reveals a strong soft component (~ 0.5 keV) significantly cooler than typically found in K- 
type stars of evolutionary stage similar to that of the HD 104237-A companion (0.7 — 1 keV). 
We interpret these findings in a possible scenario where the soft X-rays are produced mainly 
by the HAe star and the hard emission by the later-type companion. The analysis of the 
HETGS spectrum provides additional diagnostics, in particular through the He-like triplets: 
the line ratios of the Mg xi and Ne ix triplets suggest the presence of plasma at high densities 
of about 10^^ cm~^, possibly indicating accretion related X-ray production mechanism. 

The observed X-ray emission of the other sources is typical of pre-main sequence stars of 
similar ages and spectral types, except for the M-type T Tauri star HD 104237-D for which 
we find evidence of remarkably soft emission, with temperature of 3 MK reminiscent of the 
X-ray emission of the classical T Tauri star TW Hya whose emission is attributed to shocks 
in the accreting plasma. 

This research was supported by SAO contract SV3-73016 to MIT for support of the 
Chandra X-ray Center, which is operated by SAO for and on behalf of NASA under contract 
NAS8-03060. This research was based on Chandra Cycle 7 Guaranteed Time Observation 
program selected by Claude Canizares, whom we thank for providing the opportunity to 
analyze these data, and for useful discussions and comments on this manuscript. 
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Table 9. Velocity from Line Shifts Measurements 



Ion 


A [A] 


Wrad 


[km s-i]^ 


Si XIV 


6.183 


0.0 [- 


-140, +240] 


Mgxii 


8.422 


76 [- 


-61, +220] 


Mgxi 


9.169 


120 [- 


-100, +220] 


Nex 


12.135 


0.0 


[-96, +85] 


Fe xvn 


15.014 


43[- 


-52, +100] 


Fc xvn 


17.096 


41[- 


-37, +110] 


Ovni 


18.970 


43 [- 


-31, +130] 



''The 99% confidence intervals are 
listed in square brackets. 



